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Abstract 


We have studied the annihilation of galactic positrons in order to 
evaluate the probabilities of various channels of annihilation and to 
calculate the spectrum of the resulting radiation,, The narrow width 
(FWHM < 3,2 keV) of the 0,511 MeV line observed from the Galactic 
Center (Leventhal et al, 1978) implies that a large fraction of positrons 
should annihilate in a medium of temperature less than 10^K and ionization 
fraction greater than 0,05, HII regions at the Galactic Center could be 
possible sites of annihilation. 
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I« INTRODUCTION 

Leventhal, McCallum and Stang (1978) have recently reported observation 
of positron annihilation radiation from the Galactic Center using a balloon 
borne germanium detector. The observed line is at 510,7 +0,5 keV, and its 
full width at half maximum (FWHM) is less than 3,2 keV, There is also some 
evidence for the three-photon continuum from triplet positronlum annihilation. 
The 0.511 MeV line was previously seen from the solar flares of 1972, August 
4 and 7 (Chupp, Forrest and Suri 1975), but because the lower energy resolu- 
tion of the Nal detector used, only an upper limit of about 40 keV could be 
set on the line width from this observation. 

Depending on the tempera ture and density of the ambient medium, positrons 
and electrons can either annihilate directly or form positronlum. The 
importance of positronlum formation in the interstellar medium was pointed 
out by Steigman (1968), by Stecker (1969), and by Leventhal (1973), and positron 
annihilation in solar flares, both direct and via positronlum, has been 
treated by Craraiell et al, (1976). Positronlum can form either in the singlet 
state which annihilates into two 0,511 MeV photons, or in the triplet state 
which decays by 3 photon annihilation. 

Positronlum is formed by both radiative recombination with free 
electrons and charge exchange with neutral hydrogen (atomic or molecular). 

Charge exchange with heavier ions is much less important than these 
processes (Crannell et al. 1976). Once formed in a particular spin state, 
positronlum in the interstellar medium annihilates from the same state, 
because the lifetimes of both singlet and triplet positronium (10 ^ sec 
and 10 ^ sec, respectively) are much shorter than typical collision times 


with interstellar gas. 
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Sources of galactic positrons are expected to be energetic particle 
interactions , long-lived radioactive nuclei from supernova explosions, and 
pulsars (see review by Lingenfelter and Raraaty 1978). All these sources 
produce positrons at relativistic energies e 

Except for ultrarelativistic positrons, the probability for annihilation 
in flight is negligible until the positrons slow down to energies of the order 
of several hundred eV 0 (Ultrarelativistic positrons have a small probability 
<•* 10%, for annihilating in flight, but these annihilations produce only a 
broad continuum which in general cannot be observed above other continua,) 

The purpose of the present paper is to investigate the annihilation process 
at energies below a few hundred eV„ At these energies positrons either form 
positronium in flight or thermalize 8 The importance of free electrons in 
thermalizing positrons was first pointed out by Crannell et al. (1976). 

In Section II we discuss the energy loss processes of positrons in a 
partially ionized medium, and the charge exchange cross section with neutral 
hydrogen; in Section III we calculate, using a Monte Carlo simulation, the 
probability of positronium formation in flight in a partially ionized gas 
as a function of temperature, density and degree of ionization; in Section IV 
we present the rates of direct annihilation and positronium formation of 
thermal positrons as a function of the temperature; in Section V we calculate 
line shapes resulting from positronium formation in flight and after 
thermalization; in Section VI we discuss the implications of the observed 
line width on the site of positron annihilation in the Galaxy; and we 
summarize our results in Section VII* 
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11+ Energy Losses and Charge Exchange of Positrons 

When a positron is injected into a partially ionized medium, it loses 
energy to the ionized component by exciting plasma waves and to the neutral 
component by the excitation and ionization of hydrogen atoms and moleculese 
The energy loss rate in the plasma, measured in eV/cm, is given by (Book 
and Ali 1975) 

-13 , _i 

dE/dl - 1.3x10 n (M(E/kT) - M* (E/kT)) E -(nA, (1) 

e 

-3 

where E is the positron energy in eV, n fi is the electron density in cm , 

T is the temperature of the electrons in the plasma, 

M(E/kT) - 2//tt J^ kT dx x 1/2 e~ X , (2) 

o 

A * (kT/4Tm e e 2 ) 1 ^ 2 /max(2e 2 /(ra u 2 ), h/(m u)), (3) 

c C C 

u » (2E/m ) 1/2 - (8kTArm )) 1/2 (4) 

e e 

and M* is the derivative of M with respect to its argument „ 

Unlike in a plasma where the energy loss process is essentially continuous, 
in a neutral medium the positron loses a significant fraction of its energy 
whenever it excites or ionizes an atom or molecule. The cross sections for 
ionizing atomic hydrogen and exciting its 2p level from the ground state 
by electrons are shown by solid lines in Figure 1 (Kieffer 1969), The 
cross sections for exciting other levels are small in comparison to that for 
the 2p level (Kieffer 1969, Moisewitch and Smith 1968), and are not expected 
to contribute significantly to the energy loss. We assume that the cross 
sections of these processes for positrons are the same as for electrons. Also 
shown by a solid line in Figure 1 is the cross section for positronium forma- 
tion by charge exchange with atomic hydrogen (Drachman, Omidvar and McGuire 1978) 

The dashed lines in Figure 1 are positron cross sections in molecular 
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hydrogen, The total cross section is from laboratory measurements made by 

Kauppila at al„ (1977) 0 (Similar measurements by Coleman, Griffith, and 

Heyland (1974) give somewhat lower values), The ionization cross-section 

is from the compilation of Kieffer and Dunn (1966) » For the charge exchange t 

cross section above 50 eV, we use the calculation of Sural and Mukherjee 

i 

(1970), while below 20 eV, we have estimated this cross section by sub- 
tracting from the total the excitation and dissociation cross section given 
by Rescigno et al, (1976) and the elastic scattering cross section of Hara 
(1974) c We then connected the point at 20 eV with that at 50 eV as shown 
in Figure 1, We assume that above 20 eV, the difference between the total 
cross section and the sum of those for ionization and charge exchange is 
the excitation and dissociation cross section of molecular hydrogen by positronso 
In atomic hydrogen, the energy lost in an excitation collision is 10,2 eV„ 

For ionization, the mean energy of the ejected electron is 1/4 of its binding 
energy, and the distribution of the ejected electron velocities can be approx- 
imated by a Gaussian with standard deviation 0,382 otc (or ■ 1/137) (Omidvar 
1965) o In molecular hydrogen, we assume an energy loss of 12 eV in an 
excitation and dissociation collision (Herzberg 1950) and we take the same 
distribution of the ejected electrons as in atomic hydrogen except for the 
different binding energy D 
III, Positronium Formation in Flight 

We consider a positron of energy E q injected into a partially ionized 
hydrogen gas of temperature T and ionization fraction X ■ n / (n 4n ) , where 

n and n„ are the electron and neutral hydrogen densities,. If the temperature 
e H 

of the ambient medium is of order 10^K or larger, the positrons will 
thermalize before forming positronium. If, however, the temperature is lower, 
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then some fraction or the posltrona from positronium in flight by charge 
exchange vith neutral hydrogen. To evaluate this fraction, we have 
performed a Monte Carlo simulation in atomic hydrogen with ionization 
fraction X and T < lo\ e Slowing down of positrons in molecular hydrogen is 
treated separately below 0 

We then allow the positron to lose energy continuously to the ionized 
component until it makes an inelastic collision with a hydrogen atom* The 
energy, E^, at which this collision takes place is determined from 

^Tfr 1 - (5) 

E 1 

where R is a uniformly distributed random number, dE/d-t is given by equations 
(1) through (4), and a is the sum of the ionization, excitation and charge 
exchange cross sections in atomic hydrogen shown in Figure l e The probability 
that any one of these processes occurs at E^ is proportional to its cross- 
section at this energy „ If positronium Is formed, we remove the positron 
from the calculation because positronium in both the singlet and triplet 
states is expected to annihilate in a time much shorter than its collision 
time in the interstellar medium*, In an excitation or Ionization collision 
on the other hand, the positron loses energy as discussed in Section II*, 

After an excitation or ionization we repeat the process specified by equation 
(5) with a new value for E q given by the difference between and the 
energy lost in the collision, and we continue this procedure until either 
positronium is formed or the energy falls below the positronium formation 
thresholds 

The fraction of positrons which form positronium in flight before 

thermalizing, f , is shown in Figure 2 as a function of X for two cases: 
ps 

-3 13 -3 4 

n e * 0.1cm , T ■ 8000K and n e ■ 5x10 cm , T - l.lCxlO K. The parameters 
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of the first case could be representative of the war® component of the 

interstellar medium (McKee and Ostriker 1977), while those of the 

second case are essentially the same as the values considered 

by Crannell et al« (1976) for solar flares * The fraction of positrons 

forming positronium in flight increases with decreasing ionization fraction; 

for the Interstellar case and X < 0 e 05 more than half the positrons form 

positronium before they thermalisse* If n -0, about 95% of the positrons 

© 

form positronium in flight,, We also see that f Increases with increasing 
electron density because of the reduction of the energy loss rate to the 
plasma due to the reduction of the Coulomb logarithm, whose argument is given 
in equation (3)« We have calculated f for temperatures other than shown 

4 

in Figure 2, and we find that it varies by less than 20% for T < 6x10 K e For 
higher temperatures, f is expected to be small in any case because of the 
absence of neutral hydrogen. 

We have compared the results of Figure 2 with the calculations of 

Crannell et al*, (1976) and we find that for X * 0*5, our f is smaller by 

ps 

about a factor of 2 than theirs (see Figure 4 in Crannell et al e 1976)* 

This results from an erroneous factor of 2 in the Fokker-Planck equation 
that was used by them to treat the slowing down of positrons in a partially 
ionized plasma (G- Joyce, private communication 1978) 9 For X - 0*09, the 


discrepancy is much smaller because in this case, energy loss to the plasma 

is less important relative to charge exchange,, We wish to note tluit for 

low values of X, the effect of the neutrals on the energy loss becomes 

important; for example, for X = 0,09, neglecting the neutrals would Increase 

f by about 20% a 

ps 

We do not expect molecular hydrogen in a phase of the interstellar medium 
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with high temperature and high ionisation fraction. For molecular hydrogen, 

therefore, we only evaluate f for n ■ 0. Using the same Monte Carlo 

simulation as for atomic hydrogen and the cross sections of Figure 1, we 

obtain f * 0.93. This value is slightly lower than f in atomic hydrogen 
ps P 8 

at n g ■ 0 because of the higher threshold for positronium formation In the 
molecular case. We note, however, that this result is somewhat uncertain 
because of the cross sections that we have assumed. 

IV. The Fate of Thermal Positrons 

Having evaluated the fraction of positrons that form positronium in flight, 
we now consider the fate of those positrons which thermalize before annihilation. 
Thermal positrons will eventually form positronium or annihilate directly, and 
both these processes can occur with either bound or free electrons. Positronium 
formation with free electrons is radiative recombination. Radiative recom- 
bination and direct annihilation rates with free electrons were calculated 
by Crannell et al. (1976); we have extended these calculations to lower 
temperatures, and the results are shown In Figure 3. The direct annihilation 
rates with atomic electrons shown in this figure was obtained by using the 
cross section for direct annihilation, equation (3) in Crannell et al. (1976). 

We let the relative velocity between the electron and positron be ore; this 
2 2 

yields R, /n„ * 2n r c. 
aa rt o 

The rate of charge exchange of thermal positrons with residual neutral 
hydrogen can be calculated by assuming that the positrons are in thermal 
equilibrium with the an&ient medium. This is justified because even though 
annihilation and positronium formation tend to modify the positron distribution, 
the rate of reestablishing thermal equilibrium is much taster than the sum 
of the rates of all annihilation processes. Therefore, the rate of charge 


8 


exchange of thermal positrons can be written as 

*« /n » ' Cm n » <E) °« (E) v(E)dE - (6> 

where n^ is the Maxwell-Boltzmann distribution, and cr ce is the charge 
exchange cross section with atomic hydrogen shown in Figure 1, We have 
evaluated equation (6), and the results are shown in Figure 3. As can be 
seen, this rate is a very strong function of temperature, but even at its 
maximum, it Is lower than the thermal izat ion rate (Spitzer 1962). 

From Figure 3 we see that for temperatures greater than 10^K, the domi- 
nant annihilation process is direct annihilation with free electrons. At 
lower taqaeratures, radiative recombination becomes more inportant than 
direct annihilation, but depending on the amount of residual neutral 
hydrogen present, charge exchange may in fact be the most important process. 

For example, at T * lO^K, the rate of charge exchange is greater than that 

-6 

of radiative recombination if n^/n^ > 5x10 . Since for temperatures less 

t. in 10^K, we expect larger neutral hydrogen concentrations, charge exchange 

4 

should be the dominant process down to temperatures around 10 K« At 8000 K, 
R ce /n^ ■ R rr /n e , so that if, for example, the :edium is half ionized (e.g. 
the warm ionized component of the interstellar medium, McKee and Ostriker 
1977), the rates of charge exchange and radiative recombination are about 
equal. At lower temperatures, radiative recombination should remain the 
dominant process as long as there is an appreciable concentration of free 
electrons e However if becomes small, a large fraction of the positrons 

form positronium in flight, as can be seen from Figure 2- 

From Figure 3 we also see that there is a finite probability for 
direct annihilation even when positronium formation Is the dominant process. 
For example , at T ■ 8000K and *n^, about 7% of the positrons will 
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annihilate directly (mostly with free electrons), 

V , The Ener ,j * Spectrum of Positron Annihilation Radiation 

In this section, we evaluate the energy spectrum of positror* annihilation 
radiation resulting from posit ronium formation in flight and from the anni- 
^ hilatlon of thermal positrons. We first consider annihilation in flight, 

4 Using the Monte Carlo simulation described in Section III, we calculate the 

photon spectrum in molecular hydrogen with n ft ■ 0 and in atomic hydrogen 
with X ■ 0 and 0,1, As can be seen from Figure 2, for more ionization 
only a small fraction of the positrons fora positronium in flight. 

The probabilities for forming positronium in the singlet or triplet 
state are 0,25 aid 0,75, respectively. Annihilation from the singlet 
state produces two photons of 0,511 MeV in the rest frame of the positronium. 

In the frame of the galaxy, the energy of a photon from such an annihilation is 

1/2 

• e 0 (l-0cose*)/(l-0 2 ) , (7) 

where * 0,511 MeV, £ is the velocity of the positronium, and 8 * is 
the angle between the positronium velocity and the photon direction in the 
positronium rest frame. We assume that the photon direction is isotropic 
in this frame. Annihilation from the triplet state produces three photons 
whose energies can also be obtained from equation (7) provided that £5 is 

ft 

replaced by e which is chosen according to the distribution given by Ore 
and Powell (1949), 

The results are shown in panels a, b, and c of Figure 4, There is little 
difference between the line shapes for the atomic and molecular cases with 
n ■ 0„ The full width at half maximum for both these cases is about 
* 6,5 keV, The line shape for 10% ionization fraction has a flatter top 

than those for n e * 0, because in this case the rapid energy losses cause the 
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spectrum of the positronium formed to cut off sharply at low energies. 

The width of the line for X * 0,1 is about 7 keV „ The effect of the three 
photc annihilation is manliest in a slight asymmetry in the line shape and 
a low “.nergy tail. 

In the case of annihilation of thermal positrons, the shape of the 
0,511 KeV line from both singlet positronium annihilation following radiative 
recombination and direct annihilation with free electrons is expected to 
be a Gaussian of FWHM •* 0,011 HceV) (T(K)) 1/2 (Cranncil et al, 1976), For 

4 

example, if T ■ 10 K, the width is 1,1 keV, For direct annihilation 
with atomic electrons, the line shape is determined by the motion of the 
bound electrons and we expect a FWHM of about 7.5 keV, 

The photon spectra, qfe^), resulting from charge exchange cf thermal 

positrons was calculated as follows: for annihilation from the singlet state 

_ 1/2 

q(€^) - j dE >i mB (E) a ce (E)v/((E-6.8eV)(E-6.8eV + 4e 0 )) , (8) 

E 1 

where E^ 3 (&y^e o ) 2 /e^ + 68 eV, while for the triplet, 

q(e ) - r dE it_(E) o (E)v/((£-6.8eV)(E-6.8eV + P fx)/x C$) 

Y \ 2 "* ce X 1 

where 

E, - { E 1 f0r £ V >E 0 

6.8eV for e £ e 0 . 

Y ° 

For equation (8), P fc (x) is the distribution of Crre and Powell (1949) with 

ft 

x • /c q , ■ £y/(t Q (l +&)) „ mid x 2 * min (1, t;,y/(y£: o (l-3)) , where y and 
£ are the Lorentz factor and velocity (divided by c) of the positronium atom 
formed, corresponding to (E -6-8 eV) , 

By combining the results of equations (8) and (9), we have plotted in 
panel d of Figure 4 photon spectra from thermal positron annihilation for 
T * 8000K and 50000K. The FWHM's for these two cases are about 1.5 and 3.2 
keV, respectively. 
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VX» Diwigglcg 

As we have pointed out in the Introduction, the FVHM of the 0.5U MeV 
line fro* the Galactic Center was observed to be less than 3.2 keV (Leventhal 
et al. 1978) . When compared with the line shapes shown in Figure 4, this 
width is inconsistent with positronium formation in flight, and therefore 
requires that at least some of the positrons thermalize before they annihilate. 
The temperature of the medium in which they thermalize should not exceed about 
lO 5 *, since otherwise the line would be broader than observed, and the medium 
should be substantially ionized, so that a large fraction of the positrons 
thermalize (see Figure 2). 

If we consider the three component model of the interstellar medium 

(McKee and Ostriker 1977), then the most acceptable sites of annihilation 

are the warm neutral and the warm ionized media, whose temperatures are 

about 8000K and ionization fractions are 0.15 and 0.68, respectively. For 

this temperature, the width of the line is only ~ 1.5 keV, and, as can 

be seen from Figure 2, 80% and 98% of the positrons thermalize for these 

choices of ionization fractions. The bulk of the positrons should not 

annihilate in cold clouds even though most of the interstellar raa^tar is 

beleived to be in them. Since n e /n^ Is very low in clouds (< 10 Guerin 

et al. 1977), the line width would be larger than observed. Neither should 

they annihilate in the hot cong>onent whose temperature is about SxlO^K and 
-3 3 

density 3x10 cm ; however, we do not expect much annihilation in this medium 

g 

in any case, because the annihilation time is on the order of 3x10 years, 
and in this time, the positrons are expected to encounter the denser clouds 
in which their annihilation time is much shorter. 

Other possible annihilation sites are youug supernova remnants. For 
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4 

example in the filaments of the Crab Nebula, n g «* n H and T g < 10 K (Davidson 
and Tucker 1971). Supernovae are also likely sources of positrons through 
all the processes mentioned in the Introduction. 

The 0.511 MeV line was observed from the direction of the Galactic Center, 
but because of the vide opening angle of the detector (~ 15°) the field of 
vies included a large portion of the Galactic Caster region extending to 
the expanding er* of 135 km/ sec (Mezger 1974). If the positron source, 
however , is located at the nucleus of the Galaxy, the positrons would thermalize 
and annihilate in the extended BIX region in the nuclear disk, even though 
most of the diffuse matter is in the 270 pc molecular ring (Scoville 1972) 
encircling this region (Hezger 1974). The ionized region has an electron 
density of 15 to 30 cmT'* and temperature less than 10^K. From equations (1) 
through (4), we «nq>ect positrons with Initial energies of several hundred keV 
to traverse a total distance of about 500 pc in about 2000 years before they 
thermalize. Because of scattering, however, the rectilinear distance could 
be much smaller, and therefore the positrons would annihilate before reaching 
the molecular ring. The annihilation time in the ionized medium is about 

3 

10 years. For the temperature of this medium, the width of the line is 
less than about 1.5 keV, consistent with the observed upper limit of 3.2 keV. 

VII . Summary 

We have investigated the annihilation of galactic positrons, and the 

formation of the 0.511 MeV line and the accompanying 3-photon positronium 

continuum. In a cold and neutral medium, about 95% of the positrons form 

positronium atoms in flight and these atoms annihilate before undergoing 

15 -3 

further collisions if the density of the medium is less than ~ 10 cm . 

The shape of the 0.511 MeV line from positronium formed in flight is broad, 
with FWHM =- 6*5 keV (Figure 4). The width appears to be inconsistent with 
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the observed width of the 0.5U MeV line from the Galactic Center {< 3.2 keV), 
Leventhal et al. 1978). 

In a warm, partially ionized gas of temperature less than several 

4 

1 times 10 K, only a fraction of the positrons form positronium in flight 

temperatures, essentially all the positrons thermal ize before forming 

positronium or annihilating directly. 

Hie rates for the various processes leading to positron annihilation 

are shown in Figure 3. For temperatures less than ~ 10^K, positronium 

formation is the dominant annihilation channel, while for higher temperatures 

free annihilation dominates. The width of the 0.511 MeV line from thermal 

5 1/2 

positron annihilation is about 3 keV for T“10 K and varies as T . The 
observed width from the Galactic Center implies a temperature less than this 
value, and a degree of ionization larger than about 52, for which half the 
positrons thermalize. For these conditions more than 90% of the positrons 
annihilate after forming positronium, and this is consistent with the possible 
evidence for 3-photon continuum found by Leventhal et al. (1978). 

Among the various galactic sites, the warm component of the interstellar 
medium, filaments or knots in supernova remnants, and the extended HII region 
in the nuclear disk could be sites for forming the 0.511 MeV line from 
positron annihilation. 
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Figure Captions 

Figure 1* Cross sections for charge exchange, excitation and ionisation 
In atomic hydrogen (solid lines), and the total cross section, 
cross sections for ionization and charge exchange in molecular 
hydrogen (dashed lines) as a function of positron energy. 

Figure 2. The fraction of positrons forming positronium atoms before 
thermalizing, by charge exchange with atomic hydrogen as a 
function of the ionization fraction of the gas for two sets 
of parameters. 

Figure 3. The rates (per unit target density) at which thermal positrons 
form positronium by charge exchange with neutral H (R ce /ny) 
or by radiative recombination with free electrons (R rr /n e ), 
and annihilate directly with free electrons (Rd a / n e ) or with 
bound electrons as functions of the gas temperature 

Figure 4. The calculated profiles of the 0.511 MeV positron annihilation 

line. Panels (a), (b), and (c) show the profile from positronium 
formation in flight in neutral molecular, neutral atomic, and in a 
mixture of 90% atomic and 10% ionized hydrogen, respectively. 

Panel (d) shows the profile from charge exchange after thermalization 
of the positrons with residual neutral hydrogen in the 10% ionized 
case. All of these spectra include the contribution of the 
3-photon annihilation of triplet positronium. 
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